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Abstract 

We consider the scaling behavior in the critical domain of superconductors 
at zero external magnetic field. The first part of the paper is concerned 
with the Ginzburg-Landau model in the zero magnetic field Meissner phase. 
We discuss the scaling behavior of the superfluid density and we give an 
alternative proof of Josephson's relation for a charged superfluid. This proof 
is obtained as a consequence of an exact renormalization group equation for 
the photon mass. We obtain Josephson's relation directly in the form p s ~ t u , 
that is, we do not need to assume that the hyperscaling relation holds. Next, 
we give an interpretation of a recent experiment performed in thin films of 
Y Bd'iCu-zO'j-h. We argue that the measured mean field like behavior of the 
penetration depth exponent v' is possibly associated with a non-trivial critical 
behavior and we predict the exponents v = 1 and a = — 1 for the correlation 
lenght and specific heat, respectively. In the second part of the paper we 
discuss the scaling behavior in the continuum dual Ginzburg-Landau model. 
After reviewing lattice duality in the Ginzburg-Landau model, we discuss the 
continuum dual version by considering a family of scalings characterized by 
a parameter £ introduced such that itv^q ~ where m^o is the bare mass 
of the magnetic induction field. We discuss the difficulties in identifying the 
renormalized magnetic induction mass with the photon mass. We show that 
the only way to have a critical regime with v' = v ~ 2/3 is having £ ~ 4/3, 
that is, with rrih o having the scaling behavior of the renormalized photon 
mass. 
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I. INTRODUCTION 



The study of the critical behavior of the high-temperature superconductors (HTSC) 
has been initiated long ago ]1|. The discovery of this remarkable class of materials has 
opened a new perspective for the general theory of critical phenomena. In fact, the situation 
for the HTSC is much more favorable than for the conventional metallic low-temperature 
superconductors. The reason is that the size of the critical region in the HTSC is large enough 
to make its critical properties experimentally accessible using today techniques. In the tested 
temperature region it has been obtained that —0.03 < a < for the specific heat exponent, 
v ~ 0.67 for the correlation lenght exponent and the amplitude ratio A + /A~ pa 1.065 
These values are consistent with those found for 4 He. This means that the critical 
region probed corresponds to an uncharged 3D XY universality class M. In this situation 
the superfluid density, p s , scales as p s ~ A -2 , where A is the penetration depth. The 
experimental result v pa 0.67 follows from a direct measurement of A. The measured value 
of the penetration depth exponent is v' pa 0.33 and from Josephson's relation p s ~ £ -1 the 
result v pa 0.67 follows. Recent zero field experimental results obtained using very clean 
crystals of Y Ba2Cu^Ojs (YBCO) by Kamal et al. 0, confirm early measurements of the 
penetration depth near T c , giving v' = 0.33±0.01. These measurements estimate the critical 
region as being nearly 5 K wide. While this seems to be true for three-dimensional crystals, 
it does not seem to be the case for thin films. Indeed, recent measurements of A in thin 
films of Y Ba 2 Cuj,0'j-i l performed by Paget et al. || display a critical regime where the 
3D XY behavior is absent and v' = 1/2, that is, a mean field like behavior. The critical 
region reported by them is however only 0.5 wide. The same result v' = 1/2 has been 
obtained earlier for thin films by Lin et al. [[/]] . However, the result of Lin et al. followed as a 
consequence of taking in account the finite size of the sample, otherwise the result v' pa 0.33 
of Kamal et al. [||] follows. We stress that none of these experiments were performed in 
the charged critical region. The true charged transition corresponds to a very small critical 
region and is presently inaccessible to the experimental probes. From duality arguments in 
the lattice model, it is obtained that the charged phase transition should correspond to an 
"inverted" XY behavior ||. For this situation it has been argued by Herbut and Tesanovic 
that p s ~ A -1 which implies that v' = v. Recent numerical study in the lattice model 
confirms this prediction and gives v' = v pa 0.67 ~ 2/3. 

The result v pa 0.67 ^2/3 has been obtained in recent years using perturbative renormal- 



ization group (RG) methods in fixed dimension d — 3 |I 1 |14| . RG calculations performed 
in the seventies on the basis of the e-expansion leads to a flow where the charged infrared 
stable fixed point is absent if the number of components n of the order parameter is less 
than 365.9 [15]. Since the physical case corresponds to n = 2, we have that no second order 
phase transition is predicted in that calculation. The authors of ref. [|15] have predicted 
that a weak first order transition takes place. Further calculations using dimensional reg- 



ularization confirmed this scenario even up to 2-loop order [16]. This picture seems to be 



appropriate for superconductors in the type I regime. However, for the type II regime this 
result is shown to be an artifact of the e-expansion [H]. The e-expansion can be improved 
by doing a Pade-Borel resummation In this case a charged infrared stable fixed point 



is found and it is obtained that v pa 0.771 [IS]. Non-perturbative calculations on the basis 



of the Wilson RG gives v pa 0.53 by truncating the average action (the Legendre transform 
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of the Wilsonian effective action) in |0| 4 and v rs 0.58 with a truncation in |0| 8 [19]. Thus, 



the perturbative RG in fixed dimension seems to give more interesting values of v though 
it is less controlled. The non-perturbative calculations using Wilson RG of ref. [HJ are also 
performed in fixed dimension but it seems to be very difficult to make further improvements 
with respect to the different truncations. 

It is in general very difficult to find reliable approximations to study the critical domain of 
a charged transition. However, since the experimentally probed critical regime corresponds 
to a crossover near the neutral XY universality class, theoretical studies of scaling behavior 
are often performed using a Ginzburg-Landau (GL) model coupled to an external magnetic 
field H and neglecting gauge field fluctuations |[20|| . When the order parameter critical 
fluctuations are taken in account it is possible to study the different critical regimes in a 
phase diagram in the H — T plane. The situation is particularly interesting for the HTSC 
where many new physical effects like vortex-lattice melting occurs P j2l|j2^| . 

In this paper we will consider the scaling behavior in the charged critical domain in 
superconductors at the zero magnetic field Meissner phase. In the first part of the paper 
we will discuss the scaling behavior of the Ginzburg-Landau model. We use an exact RG 
equation for the photon mass to rederive the Herbut and Tesanovic result v' = v. Since 
our derivation does not use the Josephson relation P3|, we obtain as a consequence the 



Josephson relation for a charged superfluid directly in the form p s ~ t u , while in the original 
Josephson's paper this form follows only if it is assumed that hyperscaling holds since he has 
obtained actually that p s ~ t 2 ^~ w . Next, we give an alternative interpretation of a recent 
experimental result of Paget et al. obtained using YBCO thin films . 

The second part of the paper is concerned with the continuum dual version of the GL 



model [pq, p] , p9|1 . We review the lattice duality for the lattice GL model and the meaning of 
"inverted" XY behavior |J. Then we discuss the scaling behavior of the proposed continuum 
dual version of lattice duality. For this end, we consider a family of scalings where the 
parameter ( > is introduced through m| ~ ft, with rrih,o being the mass of the magnetic 
induction field. We obtain the scaling behavior for some relevant values of £ by assuming 
that the renormalized counterpart of m^o is the renormalized photon mass of the GL model 
in the Meissner phase. We show that the only way to obtain v — v' « 2/3 is using ( « 4/3. 
This means that m^ o has the same scaling behavior as the renormalized photon mass of 
the GL model. This analysis helps us to understand the difficulties in describing a correct 
"inverted" XY behavior with the continuum dual model. 



II. SCALING IN THE GL MODEL FOR T <T C 



In the following we will assume the existence of an infrared stable fixed point. This 
should be true at least in the type II regime, as lattice results have convincingly shown 
P JT0|J40|] . The existence of an infrared stable fixed point has also been established in recent 
years directly in the continuum model by using mainly RG techniques P,pT]-- [i4] , p!9| , [T7| , p!8| . 
The bare action for the d = 3 GL model in the zero external magnetic field Meissner phase 
is given by 



o,0o; m ,u ,e 



1 



d 3 x 



^(V x A ) 2 + |(V-ze o A o )0o| 



— m 
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where the subindex denotes bare quantities and wig > 0. Here wig ~ t, t being the reduced 
temperature. The renormalized action, Sr, is defined in a standard way [24|, that is, 



S R (A,frm 2 ,u,e) = S(^ /2 A, zf^ zf Z^m 2 , Z V Z;\ ^ 1/2 e), (2) 

where the quantities without the zeroes are renormalized and we have defined the corre- 
sponding renormalization constants ||25|| . The renormalization constants above are the same 



as in the Coulomb phase, as dictated by the Ward identities [EBJ. We assume that the corre- 



lation functions are evaluated in a renormalized gauge, the so called R a gauge [^6] , and that 
the Coulomb gauge limit has been taken in the end. We cannot fix the gauge a = from 
the very beginning because the unphysical fields constituting the R a gauge would become 
massless, generating infrared divergences. The unitary gauge, on the other hand, does not 
have this problem since only physical fields are present. However, it is not renormalizable, 
even in d = 3. 

In order to obtain the RG equations we must define a scaling, that is, the choice of the 
mass scale controling the flow while specifying the variables which are kept fixed in this 
process. Here we note that we have two dimensionful couplings, namely, e 2 and u, both 
having dimension of mass. In the Meissner phase the photon acquires a mass, m^, which 
together with the mass, m, define two mass scales in the problem. An useful dimensionless 
parameter is the ratio between these two masses, the Ginzburg parameter, k = m/mA = 
(life 2 ) 1 / 2 . Since the critical point corresponds to m = 0, the renormalized mass m is a 
good mass scale to control the flow. Thus, we will choose it as the fundamental mass 
scale. Let us define in this way the dimensionless couplings u = um d ~ A = Z^Z~ 1 uom d ~ A 
and e 2 = e 2 m d ~ A = Z \e 2 {m d ~ A . Note that although the GL model has been written in 
fixed dimension d = 3, we have defined the dimensionless couplings in arbitrary dimension 
2 < d < 4. In order to complete the definition of our scaling we must specify the parameters 
that will be kept fixed. The standard choice corresponds to differentiate the renormalized 
dimensionless couplings with respect to ln(m) keeping the bare couplings and the ultraviolet 
cutoff fixed. Note that in this process neither m nor m^o are kept fixed. We define the RG 
functions 

d\nZ A 

Va = m^— (3) 



d\nZ^ 



(2) dhlZ? 

% = m ^r (5) 

At the infrared stable fixed point T)a and rj^ gives respectively the anomalous dimensions of 
A and <\>. The anomalous dimension of |0| 2 is given by the fixed point value of m — rj^. 
The fixed point values e 2 and are determined from the equations (3^ = mde 2 /dm = 
and fin = mdu/dm = and we take the solution corresponding to the infrared stable fixed 
point. The beta funcion (3^ is given exactly by 
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= (vA + d- 4)e 2 . 



(6) 



An immediate consequence of the above equation is that a charged fixed point corresponds 
to r/\ = 4 — d ||^ |rByiiyr5| ; V% being the fixed point value of tja- From Eq. (^j) and the 



definition of k we have the following exact equations: 



9k 2 
dm 



K 



u 



+ A-d-r] A 



(7) 



m _ = mA ^-2 + ^- T J. (8) 

From Eq. (§) we obtain easily that near the phase transition (that is, near the charged 
infrared stable fixed point) the photon mass scales as itla ~ m. Since m = £ _1 and = 
A" 1 , A being the penetration depth, we obtain A ~ £ implying v' = v. This is a rederivation 
of the result of Herbut and Tesanovic ||. Note that in our derivation no use has been made 
of the Josephson relation. Since m\ = e 2 p s and from Eq. (|6|) e 2 ~ m 4 ~ d = £ d ~ 4 , we obtain 
that p s ~ £ 2 ~ d ~ fid-- 2 ) This constitutes a renormalization group proof of the Josephson's 
relation for the charged superfluid. In the original Josephson's paper this relation is proved 
for an uncharged superfluid and given in the form p s ~ t 2f3 ~ vu (77 is the fixed point value 
of 770). This form follows easily by noting that p s =< |0| 2 >= Z^ 1 < |0 O | 2 >. Since near 
the critical point ~ m' n ~ t vr) and defining j3 through < \<f)o\ 2 >~ t 2 @, it follows that 
p s ~ t 2l3 ~ ur] . Thus, in this last argument leading to Josephson's relation it does not matter 
if the superfluid is charged or not. Our derivation made directly for the superconductor 
implies therefore 

2/3-771/ = v(d- 2). (9) 



As pointed out by Fisher et al. [£8[] in the context of uncharged superfluids, the relation 
with the exponent v(d — 2) instead of 2/3 — r/v holds only if hyperscaling holds, that is, 
dv — 2 — ol. Our argument shows that hyperscaling holds for a superconductor. For an 
uncharged superfluid the hyperscaling relation can be proved by using scaling RG arguments 
of the same type we used here |24[]. Of course, this type of proof is not rigorous in the 
sense that it assumes that the continuum limit of the lattice statistical mechanical model 
exists. Only the inequality dv > 2 — a can be rigorously proved [57]. A well known case 



where hyperscaling fails is mean field theory where v — 1/2 and a = independent of the 
dimension. In this case hyperscaling is fulfilled only at d = 4. 

The Josephson's relation is used experimentally to determine the value of v. Today it 
is possible to perform very accurate direct measurements of the penetration depth. The 
critical region probed is such that the gauge field fluctuations are unimportant and the 
critical fluctuations are those of the order parameter field and this means that 77^4 = 0. 
Then, from Eq. (Q), we obtain that near the phase transition m 2 A ~ m d ~ 2 ~ ^( c( - 2 ) ) that is, 
p s ~ A~ 2 ~ £"(a- 2 )_ Thus, in such critical regime 

✓ - ^ do) 
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Experiments performed in YBCO crystals give that v' = 0.33±0.01 [0,0- Using Eq. ( |10"D for 
d = 3 we obtain v m 2/3, consistent with the 3D XY behavior, that is, a 4 He like behavior. 

The situation seems to be however different for YBCO thin films. A recent measurement 
by Paget et al. || performed in YBCO thin films gives v' = 1/2, that is, a mean field like 
behavior. This result has been obtained to within 0.2 — 0.5 K of T c , indicating in this way a 
much smaller critical region as compared to the bulk YBCO However, this mean- field 

like behavior could be interpreted as a non-mean field behavior in the following way. If we 
insist in using (pi]) to evaluate v we obtain 

(11) 



d-T 

and, assuming that hyperscaling holds, 

d — A 

The above exponents are not classical and we have v = 1 and a = — 1 for d = 3. Note the 
similarity between the present critical regime with the 0(n) model for large n and 2 < d < 4. 
The above exponents are just the exact exponents for the 0(n) model at large n EM. The 



exponent v as given in ((TTJ) is also obtained in a 0(n) non-linear a-model in d = 2 + e 
dimensions and n > 2 Thus, it is possible to interpret the experiments of Paget et 

al. as correponding to a non-classical situation characteristic of superconducting thin films. 
There are, however, some possible handicaps in this point of view. For instance, it is not 
expected a so expressive change in a for thin films relative to the bulk material. Thus, if a 
remains close to zero, the critical behavior probed in ref. || corresponds in fact to mean- 
field. Another important point is the possibility of dimensional crossover behavior due to 
finite size effects arising from the thikness of the film. This situation has been extensively 
studied in superfluid 4 He films [p9|-|3T|. For superconducting YBCO thin films the situation 
is unclear because in the optimally doped case the coupling between the CuO planes is 
strong and the 2D fluctuations are probably dominated by the 3D fluctuations even for a 
small film thikness. However, in the underdoped cuprates the coupling between the CuO 
planes is much weaker and we can expect a strong influence of 2D fluctuations for sufficiently 
small film thikness. 

To conclude this section, let us discuss the scaling behavior of the order parameter. 
This is a controversial matter both from the theoretical and experimental point of view. 
The theoretical controversy has its origin in the scaling behavior of the correlation function 
Wq 2 '(x, y) =< o (x)0o(y) > a ^ l ar S e distances, |x — y| — > oo. For |p| — > and m = 
its Fourier transform Wq 2 ^ behaves as W (p) ~ IpI* 7-2 - Most calculations gives a value of 
7] in the range —1 < r\ < P,|l3,19,15-|i~8||. This does not contradicts the scaling relation 



(3 = + T])/2 (for d = 3) because — 1 < rj < imply (3 > as it should be (note that 
v must be positive). However, a negative value of r\ is pathological in many respects. For 



instance, it has been pointed out by Kiometzis and Schakel [|32) that rj < violates the 
positivity of the spectral weight in the Kallen-Lehmann spectral representation of Wq |33 



In fact, this representation of W implies that < < 1. Since near the critical point 
~ m v , we have necessarily that rj > 0. Moreover, from the |p| — ► behavior of wfj 2 \p) 
we see that rj < makes the low momentum behavior worse than before renormalization |34| . 
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This is in contradiction with the infrared stability of the fixed point. Another important 



point is that i] is the fixed point value of rj^, a quantity which is gauge dependent |L4 
Thus, we can question the physical meaning of the r\ exponent and, as a consequence, of 
the order parameter itself. This situation is unconfortable because in the Gorkov derivation 
of the GL model from the BCS theory P^| , (p is defined as being proportional to the gap 
function. In the microscopic theory the gap function has a precise physical meaning since 
it is responsible for the generation of a gap in the spectrum. The expectation value of <fi is 
thus proportional to the gap near the critical temperature. The gauge dependence of the 
superconducting order parameter has been discussed recently in \TA\ and it has been shown 
that drj^/da — > as the critical point is approached (a is the gauge fixing parameter), the 
only gauge contribution left to r\ corresponding to the Coulomb gauge, a = 0. This result 
agrees with an early analysis by Kennedy and King |36|], although these authors used a 



different, but closely related, definition of order parameter. Thus, the gauge dependence is 
actually not an issue and the point to be solved in the approximations is the negativeness 
of Tj. 



III. SCALING AND DUALITY 

In this section we will study the scaling behavior of superconductors in a continuum dual 
Ginzburg-Landau (dGL) model. The dGL model has been proposed using plausible argu- 
ments on the dynamics of a vortex gas [Tl|j3^j3^| and is assumed to be the continuum version 
of the dual GL model in the lattice 



Lattice duality studies in abelian gauge models 
P| |?7|Pm helped condensed matter theorists to obtain important conclusions concerning the 
superconducting phase transition. In particular, it has been used to establish that a second 
order phase transition should take place at least in the type II regime H. 



A. Duality in the lattice GL model 



For the sake of clarity, we set up here the arguments given in several papers 
A lattice version of the GL model has a partition function given by 



Z(/3,e) 



-r-r dOi 

A ^ 2tt 



]JdA 



ifj. 



1,1 1 



exp 



P cos(A^ - eA ifl ) - - £(A x A; 



1 1 1 



, (13) 



where A^Oi = — 9i, that is, is a lattice derivative. Note that in the above partition 
function the integral over A ifJi is over the interval (— oo, oo), meaning that the gauge group 
is 1R, a non-compact gauge group. 



In the Villain approximation Q we can rewrite fll3|) as 



Z(p,e) 





-n- dOj 


POO 


Y[dA ilx 








1 — 7T 


_ % 


J —oo 







■| J2( A A - eA in ~ 2vrm iM ) 2 - ~ XX A x Aj 

ifi i 

(14) 
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Let us introduce an auxiliary field 6j„ such that 



cxp 



--(A^ - eA^ - 2-nm^) 



/oo 
db ifl exp 
-oo 



(15) 



In the following we will neglect factors of proportionality which are generally smooth func- 
tions of (3. Thus, many equations will be assumed up to proportionality factors. By per- 
forming the 9 integrals we have 



Z(p,e) 





/•OO 






J —oo 


Y[db ilM 


i,fi 







5(A-bi)5>xp 



hi 



(1 



Applying the Poisson formula we obtain 



1,1 1 



£ ^A-n l5 exp 



• (17) 



The constraint A • rij = in ( |TTD implies that the links iii^ form closed loops. These closed 
loop are interpreted as magnetic vortices |3l],|37]]. It can be shown that the XY model 
partition function in the Villain approximation (Eq. (13) with e = 0) can be cast in the 
form 



Zxy(P) 



]~| da 



1 1 1 



^a-m^o exp 

Mia 



]T - — (A x a^) 2 + 2m a^Af, 



in 



(18) 



and we have Z(oo,e) = Zxy(e 2 /(27r) 2 ). j3 < oo corresponds to add a chemical potential 
for the loop variables in the dual XY model. In this sense we can think of the lattice GL 
model in (|14j) as a generalized dual XY model. Precisely, we have Zxy(P) = Z XY (/3,0) 
where Z' XY (/3, K) is the partition function of a generalized dual XY model given by 



Z' XY ((3,K) 



Y[da ifl 



Miu, 



V (A x a. 



+ 2m Oi/tMift 



i/j. 



(19) 



It is possible to study the phase transition of the above dual model by looking the phase 
diagram in the j3 — K plane ||. In this phase diagram the point (/5 C , 0) corresponds to the 
XY critical point, f3 c being the inverse critical temperature. From Eq. ([TJ]) with e = we 
obtain after integration of the 9 variables 



Zxy{0) = Z(P,0) = E 5 A.m l5 oexp ( ~7^E m 



(20) 



IfJL 
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On the other hand, integration of the in Eq. (jl9|) yields 



z 'xy{P, k ) = Y1 exp 



M„ 



■i 1.1 



(21) 



where the lattice Green function G behaves like |xj — Xj| _1 for |xj — Xj| — > oo. From Eqs. 
(p0|) and (|2T| ) we have that the point (0, l/2/3 c ) in the (3 — K plane corresponds to the so 
called "inverted" XY transition. By doing the integration in (jl7|) we obtain 



.2 j 



^,e) = ZW I — , — 



(22) 



which means that the lattice GL model will undergo an "inverted" XY transition for < 

e 2 <e 2 c = 4ir 2 (3 c @. 

Let us work out further Eq. fll7|) . Since A ■ rij = we can put rij = A x lj. By using the 
Poisson formula to introduce a continuum field and integrating out the A ifl y s, we obtain 



Y[dh itl 



ifj, 



S ±-™ifi exp I 



2/T tJ 2 1 



+ 2m m i^ip, 



ii-t 



(23) 



Note that putting e = in ( P3|) we obtain a partition function identical to fll8|), as it should 
be. Reintroducing the variable 6 through 



(24) 



and using the Poisson formula to convert the integral over hia in a sum over nj„, we obtain 



Z(P,e) 



"I r 

L / ^"j 



-^ £(A^ - 2ttti*) 2 - f £( A x n *) 2 
57r i/x / » 



(25) 



The partition function fl25|) is the dual representation of (|T4|) used in the numerical simu- 
lations in ref. [BJ. Note that the dual representation (Effi) has gauge group Z. We observe 



that the Z given in (|T4] ) with e = should be the same, up to proportionality constants, 
as the Z given in (25) with (3 — provided we put /3 = e 2 /47r 2 in (14). It is clear that the 
temperature plays the role of a dual charge satisfying the Dirac condition ef3~ l = 2tt. 
By performing the /i^ integration in (^) we obtain an equation analogous to (|2lD with 
= 0. The difference is that the Green function G is replaced by a massive Green function 
<5(|xj — Xj|) which behaves like e _/3e ' Xi_Xj '/|xj — Xj| for |xj — Xj| large. Thus, we can write 



- 2 1 



Z((3, e) = Z' XY f — , — J = ^A. mi ,o exp 



-2n 2 p m if iG(\yii - x^Dm^ 



(26) 



Therefore, adding a chemical potential to the loop variables in the dual XY model is equiv- 
alent to replace the massless Green function G by a massive one in a dual XY model with 
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K = 0. In this duality map we have represented the partition function of the lattice GL 
model in such a way that the loop variables have zero chemical potential. This amounts 
in replacing the massless Green function by a massive one. This means that if we consider 
a scaling (continuum) limit of the dual model, it is necessary to consider a massive gauge 
field h(x). This gauge field should satisfy the constraint V ■ h = and should be coupled 
minimally to a field ip such that \tp\ 2 represents the density of magnetic vortices. 



B. The continuum dual model 



The continuum version of lattice duality, the dual Ginzburg-Landau model (dGL) has 
been proposed on the basis of plausible arguments concerning the dynamics of a vortex 
gas [n],[8yj!J. Attempts have been made to justify it as the continuum limit of the lattice 
dual model p9|JT^| but it does not exist to date a rigorous mathematical construction of the 
continuum limit of the lattice dual model. For instance, a possible way towards motivating 
the dGL model from the lattice dual model has been proposed by Herbut |T2|]. Instead 
combining (531) and (p|) to obtain (B3), we can follow Peskin [37[ and insert in (E3I) 



1 = limexp 

t->o 




(27) 



By using (|J) and the identity 



E e 5 

m=— oo 



m +ixm 




^(x-2-kM) 2 



(28) 



we obtain 

Z(/3,e) = 



lim 



2tt_\ n/2 
t ) 



[{dh ilM 



I ft 



L( A x h,) 2 - -h 2 

2/T ; 2 1 



n 



d9i 
2tt 



E ex P 



J2Y f ( A A- 2nk ^- 2nM ^ 2 

ill 

(29) 



The limit t — generates delta functions in the integrand which, when replaced by the 
integral representation and applying the Poisson formula allows us to recover (^). It has 
been proposed in [12|] that by leaving t finite Eq. ( p9f ) is analogous to (0), with the difference 
that in ( p9[) the lattice gauge field is massive. Moreover, in virtue of (^) t can be interpreted 
as a chemical potential for the loop variables. Thus, by keeping t small but finite we can 
write the following continuum limit of (|29| ) [|11|,|T^,^,^] in terms of bare quantities: 



S 



1 . <2 %,0.2 



2 (Vxh^ + 



^0, 



(30) 



In (|30|), the bare dual charge, go, is related to the charge eo by the Dirac condition q eo = 2n. 
The field ip is usually called a disorder parameter, as opposed to the order parameter in the 
GL model. We can of course criticize the motivation of (BUI) from (f2ll) because it is based 
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on the hypothesis of a finite, though small t in (|29|). In fact, ( p9| ) represents a lattice GL 
model only if t — > 0. A small nonzero t constitutes an approximation and the continuum 
limit above should be regarded as an approximate continuum dual model. This should 
be contrasted with the exact duality map we have obtained in the lattice (up to a Villain 
approximation). The dGL model given by (|30| ) has been also motivated by arguing directly 
in the continuum but using the London limit of the GL model, that is, by assuming that 
0(x) = <j)e l6 with (f) constant [R^|. All of this is approximate and we cannot really say that 
the dual map in the continuum is exact. The construction of continuum limits is a highly 
non-trivial matter, even for a simple scalar model P5fl . 

Although the dGL model gives a respectable value for the correlation lenght exponent 



TT|,|T2|1 , it does not give v' = v when mh,o is identified with the bare photon mass rriA,o of 



the GL model | IT| . Depending on the choice of scaling (to be precised later), we can find 



a different value of z/, as will be shown in the next paragraphs. We will consider a family 
of scalings where ( > is a parameter. For each ( we find a different phase transition 
regime that must correspond to a possible critical behavior in the GL model. By considering 
this family of scalings, we will show that it is not possible to describe an "inverted" XY 
behavior with mh,o — ^Afl and, at the same time, rrih = rriA, uih being the renormalized 
counterpart of m h $. 

We define the renormalized fields h = Z h 1 ^ 2 h and ip = Z^^tpQ. From the Ward 
identities we obtain that the term m\ \\^/2 does not renormalize, implying in this way 
that m\ = Zhm\ Q . The renormalization of the remaining parameters follows easily: fi 2 = 

(z1 2 ^) -1 Z^/Zq, v = Z~ l Z 2 v G and q = q . We note that the dual charge remains bare. Let 
us define the quantity g 2 = Tn\ Q q 2 which renormalizes as g 2 = Z^g 2 . The dimensionless 
couplings relevant to the problem are g 2 = g 2 / '/i and v = v/fi. We will perform our RG 
analysis by introducing a scaling hypothesis for m^o- We assume that m\ Q oc where 
( > 0. This scaling hypothesis will define under a RG a family of scalings £f . It is useful to 
define the ratio = fi/rrih, analogous to the Ginzburg parameter k in the GL model. For 
arbitrary (>0we obtain the flow equations: 

a 2 

^ = [2-%-C(2 + 4 2) -^)]4 (31) 

^ = k + C(2 + ^ 2) -^)K, (32) 

H^ = [Vh-l + ((2 + V$ ) -V«,)}9 2 , (33) 

(2) 

where the RG functions rfl , rj^ and r)h are defined by 

(2) 9 ^ Z T ( o A , 

=»~d^' ( 34 ) 
W = A 4 qJ^ > (35) 
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dhiZ, 



(36) 



Near the phase transition // = ~ t v . As ( varies, the above equations describe 
different regimes of phase transitions in the case where an infrared stable fixed point exists. 
Let us analyse some relevant cases. 



1. Sq scaling 



For C = we have that the infrared stable fixed point for Eq. (33) is a non- vanishing 
gl satisfying rf h = r]h(g 2 ,v*) = 1, Pv(g 2 ,v*) = with {>* being the fixed point value of the 
coupling v. From Eq. (^) we obtain that near the phase transition ~ /i 1 / 2 . By defining 
the exponent of rrih through rrih ~ t Uh , we obtain that Vh = v /2. This suggests that k,i — > 
as approaching the critical point. This is the case indeed, since from Eq. ( pTj ) we obtain 
that near the infrared stable fixed point ~ fi 1 ^ 2 . From this behavior of Kd we obtain that 
the fixed point value corresponds to the value in a XY model and v = vxy ~ 2/3. The 
reason for this behavior comes from the fact that every power of the dual charge in f3% is 
multiplied by a function of which tends to zero as the critical point is approached. The 
same type of behavior has been already encountered in the litterature |i2| , |T4|| . If we identify 
rrih with the renormalized photon mass itla we obtain v' = vxy/2, the half of the value 



expected for the superconducting transition in the "inverted" XY universality class [flO 
However, this result corresponds to the universality class of the crossover regime governed 
by the neutral XY fixed point M , which is the only one we have experimental access [0.RJ . 



2. Si scaling 



We consider now the scaling Tii corresponding to £ = 1. In this case we obtain that the 
infrared stable fixed point corresponds to g 2 = 0. The present situation matchs with the 
one encountered in ]l ]J where m^o is identified with the bare photon mass m^.o- We have 
still that v = vxy- However, if = itla we obtain that Vh — v' = 1/2 and the penetration 
depth exponent is classical. As we have already discussed in section II, the value v' = 1/2 
has been measured by Lin et al. and in a recent paper by Paget et al. 0. 



3. S^ scaling with £ ~ 4/3 

Clearly the scalings S and Ex of the dGL model are not in the "inverted" XY univer- 
sality class but belong to different crossover regimes in the GL model. In order to have an 
"inverted" XY behavior we need to obtain = v' = v = uxy, after identifying with 
the renormalized photon mass in the GL model. It is easy to check that this is in fact the 
case for ( w 4/3, but in this case ro^o clearly does not corresponds to the bare photon mass. 
Rather, rrihfl behaves like the renormalized photon mass. For ( w 4/3, the infrared stable 
fixed point corresponds to a neutral dual charge just as in the scaling Si and from Eq. ( |32|) 
we obtain rrih ~ fi near the fixed point. Note that in contrast with the other two cases, now 
we have that approaches a non-vanishing fixed point value in the infrared. 
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IV. CONCLUDING REMARKS 



We have discussed in this paper some relevant aspects of the scaling behavior of super- 
conductors in a zero magnetic field Meissner phase. We have obtained in the first part of 
the paper the scaling relation v' = v as a consequence of an exact flow equation for the pho- 
ton mass. The interesting point of our derivation is that since no use has been done of the 
Josephson relation, this relation is obtained directly in the form p s ~ t u . We argued that this 
implies the validity of the hyperscaling relation for a charged superfluid since the Josephson 
relation can be derived independently in the form p s ~ t 2l3 ~ ur) . In this part some emphasis 
has been given to the relation between theory and experiment, particularly in what concerns 
the usefulness of the Josephson relation in determining the exponent v. While this is not 
new, it helped us to propose an alternative interpretation of a recent experiment performed 
by Paget et al. |J in YBCO thin films. We proposed that the correlation lenght and specific 
heat exponents are possibly non-classical and given respectively by v — 1 and a = — 1. In 
order to confirm this it is necessary to measure a directly in thin films, which is a very 
difficult task. We are aware however of the little probability of changing so dramatically 
the value of a with respect to the bulk sample value. Anyway, such a possibility cannot 
be completely discarded and we hope that some experimental effort could be made in this 
sense. 

In the second part of the paper we tried to elucidate the scaling behavior of the continuum 
dual Ginzburg-Landau model. This has been done by the introduction of a scaling hypothesis 
on the bare magnetic induction mass. Our arguments in this part of the paper show that 
the continuum dual model still deserves more reflexion. We think that the dual map in the 
continuum is not as complete as it is in the lattice. We emphasize that there is no rigorous 
argument in favor of the usual continuum model (|30|). A future theoretical perspective is to 
study different scalings in another possible choice for a continuum dual model. It consists of 
a continuum XY model where a perturbatively non-renormalizable interaction of the form 
{ip^Vip — ipVip^) 2 * s a dded. It is possible that such an interaction can be non- 

perturbatively renormalizable, like the situation encountered in the Gross-Neveu model in 
d = 3 jKJ where the non-perturbative renormalizability has been rigourously demonstrated 
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